Journal of the Physical Society of Japan 



LETTERS 



o 

(N 
C 



c3 



C 

o 

(N 
> 
On 



(N 



Quantum Thermal Hall Effect in a Time-Reversal-Symmetry-Broken 
Topological Superconductor in Two Dimensions : Approach From Bulk 

Calculations 

Hiroaki Sumiyoshi * and Satoshi Fujimoto 
Department of Physics, Kyoto University, Kyoto 606-8502, Japan 

We discuss thermal transport of two-dimensional topological superconductors (TSCs) with 
broken time reversal symmetry, which are described by Bogoliubov-de Gennes (BdG) Hamil- 



tonians. From the calculations of bulk quantities only, without refereeing to Majorana edge 
states, we show that the thermal Hall conductivity of two-dimensional TSCs in the low- 
temperature limit is quantized in multiples of which is exactly one half of the value 
^ of quantization in the case of the integer quantum Hall effect, and that this exact half- 

<D ' quantization is caused by the structure of the Nambu spinor and the particle-hole symmetry, 



which BdG Hamiltonians generally have. In the case of spinless chiral p-wave supercon- 
ductors, this result is in perfect agreement with the argument based on the Ising conformal 
field theory with the central charge c = 1/2, which is an effective low-energy theory of the 
Majorana edge states. 
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Introduction — Topological insulators (TIs) and superconductors (TSCs) are new quan- 
turn phases which are characterized by bulk topological invariants. It has been found that 



some of these invariants directly characterize the transport properties of the system. For ex- 
ample, the integer quantum Hall effect (IQHE) state is characterized by the topological num- 
ber, the TKNN (or first Chern) number C\, and this invariant appears in the Hall conductivity 
crxy, in the low-temperature limit: cr^, = C\e 2 /h, which has been derived from the linear re- 
sponse theory.^ A similar idea is also applicable to topological superconductors without time 
reversal symmetry in two dimensions (i.e. class C and D in the topological periodic tableP), 
for which topological nontriviality is characterized by the nonzero TKNN number. However, 
in this case, it is not charge transport but thermal transport that characterizes the topological 
feature, since charge is not conserved, while energy is still a conserved quantity in supercon- 
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ducting states. For instance, in the case of spinless chiral p-wave superconductors, which have 
the TKNN number C\ = ilj^Read and Green showed that by using the effective low-energy 
theory for Majorana edge states, the thermal Hall coefficient in the low-temperature limit is 
precisely given by (in the unit of k& = 1, % = 1), where c = 1/2 is the central charge 
of the Ising conformal field theory (CFT) which describes the Majorana edge stated Note 
that this value is one half of the value of IQHE states, which have edge modes described by 
the chiral Luttinger liquid theory, i.e. CFT with the central charge c = l. 5 The Read-Green's 
argument is based on the CFT for the edge states, and the transport coefficients are not cal- 
culated from the bulk Hamiltonian. This means that the bulk-edge correspondence is ensured 
only by the existence of the underlying effective low-energy topological field theory, i.e. the 
gravitational Chern-Simons theory for thermal responses. 5 However, the CFT for the edge 
state can not predict the sign of the Hall conductivity, though its magnitude up to the sign is 
correctly determined. For instance, the above-mentioned formula is derived by assuming 
that its sign is positive. Actually, the sign of the Hall conductivity depends on the bulk band 
structure. Also, it is highly challenging to derive the gravitational Chem-Simons theory from 
microscopic Hamiltonians of solid state systems. Thus, it is desirable to obtain the quantum 
thermal Hall conductivity directly from the bulk microscopic Hamiltonian without referring 
to the edge theory. In fact, the calculation of the quantized thermal Hall conductivity from the 
bulk Hamiltonian is highly nontrivial even for the case of the Chern insulator (or the IQHE 
state), as was noticed by Qin, Niu, and Shi recently.^ They pointed out that to obtain quan- 
tized thermal Hall conductivity which is associated with the TKNN number, one needs to 
extract the contributions of the energy magnetization current, since the magnetization current 
does not participate into the transport current induced by applied temperature gradient. Their 
procedure successfully gives the correct expression of the thermal Hall conductivity for the 
Chern insulator = at sufficiently low temperatures. Also this result is in agreement 
with the edge theory of the Chern insulator, i.e. CFT with the central charge c = 1 : if there 
are n-channel of the edge modes, the edge theory gives = c^n which coincides with the 
bulk result up to the sign, implying |Ci| = cn = n. 

Recently, Nomura et al. clarified the relation between the bulk TKNN number and the 
half-quantized thermal Hall conductivity in the case of two-dimensional (2D) Majorana 
fermions on the surface of three-dimensional TSCs, which is related to 2D spinless p + ip 
superconductors. 6 8 However, for more general cases of 2D TSCs such as chiral rf-wave super- 
conductors and 5-wave superconductors with the Rashba spin-orbit interaction (SOI) ElEnEH 
it has not been well understood how we can associate their TKNN number, which is an in- 
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teger, with the CFT prediction mentioned above; i.e. the thermal Hall conductivity is not 
integer-multiple of ^f, but half-quantized. To solve this issue, in this letter, we microscopi- 
cally derive the thermal Hall conductivity from the bulk Hamiltonian, without using the edge 
theory. There are two important key points for the microscopic calculation of the thermal Hall 
conductivity: (i) we should take into account the energy conservation law properly, and (ii) 
the contributions from energy magnetization must be extracted to obtain transport currents 
as pointed out by Qin, Niu, and Shi.^For this purpose, we, first, clarify the symmetric prop- 
erties which Bogoliubov-de Gennes (BdG) Hamiltonians generally possess. Next, we define 
the energy current operator which preserves the continuity equation. Finally, we calculate the 
thermal Hall coefficients of TSCs by introducing gravitational fields which act as mechanical 
forces inducing heat currents, and by using the Kubo formula with corrections from energy 
magnetization contributions. 7 It is shown that the thermal Hall conductivity is quantized in 
the low-temperature limit exactly to one half of the value in the case of Chern insulators and 
IQHE states. 

Model, symmetry, and preliminaries — In this section, we introduce a model which de- 
scribes 2D topological superconductors with broken time-reversal symmetry. The argument 
developed in this letter is applicable to large classes of 2D topological superconductors such 
as chiral 77-wave superconductors, 5 chiral d- wave superconductors, and s-wave superconduc- 
tors with the Rashba SOI under magnetic fields PCSEH 

The Hamiltonian we start with is given by 

ik = j d d 4\r)-H BdG (r,d r )$(r), (1) 

where each operator has indices (a, s). a(= +1) is an index of the Nambu space and s is a 
spin index, i^-ij(r) and \j/u(f) are defined as c s (r) and cj(r), respectively, where c s (r) (cj(r)) is 
an annihilation (creation) operator of an electron, and the BdG Hamiltionian H BdG generally 
takes the formJEEl 



HBdGasa's'{ r ,d r ) 



a' = -1 a' = +1 
a = -l K s .Ar,d r ) -A s Ar,d r ) , (2) 
a = +l A;/r,3 r ) -K* ss ,(r,d r ) , 
where K ss >(r, d r ) and A ss >(r, d r ) are operators including differential operators and functions of 
r. In the case where K(r,d r ) = -^V 2 + V(r) - p and A(r,d r ) = -A(-id x - d y ) [V(r) is a 
potential energy and A is a complex constant], the model describes spinless chiral 77-wave 
superconductors. 5 On the other hand, in the case where K{r,d r ) SS ' = {(-^V 2 + V(r))6 ss > - 
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hcr z + (T ■ g(d r )} ss > and A ss t(r, d r ) = (-iAo- y ) ss > [h is a real constant, A is a complex constant, 
cr = (cr x , cr y , <t z ) are Pauli matrices, and g(d r ) = 2A(-id y , id x , 0)], the model describes s-wave 
superconductor with the Rashba spin-orbit interaction (SOI)P^El 

For the derivation of the formula for the thermal Hall conductivity, it is useful to clarify 
the symmetry of our system. The BdG Hamiltonian is Hermitian and preserves the particle- 
hole symmetry (PHS): H BdGasa ^> = -H* BdG _ as _ a , s ,. Consequently, if f as {r) is the solution for 
the eigenvalue E, fl as (r) is the solution for the eigenvalue -E. As a result, if the system is 
periodic, the eigenfunctions for H b< ig can be written as the form f knas (r) = e' kr u knas (r), where 
Uknas(r) has the same periodicity as the system, and the eigenenergys and eigenfunctions can 
preserve the symmetries 

E kn = -E- k - n , fknasir) = fl k _ n _ as (r). (3) 

We take E nk (n = 1, 2, 3 • • • ) positive and E nk (n = -1, -2, -3 • • • ) negative, and normalize the 
eigenfuctions so that ^ as f d d rf* nas {r)f knas {r) = 1. 

Now, we obtain the expressions for correlation functions of creation and annihi- 
lation operators of Bogoliubov quasiparticles, which are utilized in the following cal- 
culations. We introduce annihilation (creation) operators of Bogoliubov quasiparticles: 
fykn = £ aj j ' d d rf k \ as {r)ij/ as {r). By using the symmetry eq. ©, we can derive the sym- 
metry of these operators ij/ kn = $_ k _ n and the anticommutation relations: ft (J ,} = 
8 k k>8nn>, {fik„Ak>n>} = {fy\ n , ft „/( = 5 k - k >5 n - n >. (Note that the latter is nonzero.) Since 
the Hamiltonian can be transformed into jk = \Y.] m E kn \^ k Ji kn = Yukn,n>^E kn ^\jt kn + 
(c - number), we get the expressions for correlation functions: 

(ftknftk'n') = f{E k 'n')8 n -n'8 k ~ k ' = (1 - f (E kn ))8 n - n > 8 k - k > , 

{^IJ/k'n') = f(E k 'n')8 im '8kk', 

(fykn$\, n ) = (1 ~ f(E k , n ,))S nn >S kk >, 

(ftnft'n') = - fi E kfn'))S n -n'Sk-kf = f(E kn )5 n - n >5 k - k >, 

where the function f(E) = l/(e^ E + 1) is the Fermi distribution function at the temperature 
T = 1//3 and the chemical potential [i = 0, and (X) is the the statical average of the operator 
X: {X) = Tr^e - ^^] /Trfe - ^]. Applying Wick's theorem, we get the expression for the 
four-point correlation function: 

+ f{E x ){\ - f{E 2 )){8 X2 ,8 2V - 8 X - V 6 2 - 2 ,), (4) 
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where Eu 6u and <Si_ 2 are the abbreviations for E kini , 6^k 2 5 nin2 and S^-^S^-^ respectively. 
Note that the last term of ©, -6i-y6 2 -2' , does not appear in the four-point correlation function 
of the systems of ordinary fermions, and is inherent for the Nambu spinor. This term plays an 
important role for the derivation of the half-quantized thermal Hall conductivity as shown in 
the following. 

Energy current operator — In general, for the microscopic argument on transport phe- 
nomena, it is important to define correctly the "charge" and "current" field operators which 
satisfy the continuity equation. Therefore, in this section, to consider thermal transport, we 
define the energy density operator and the energy current operator. To apply the linear re- 
sponse theory, we introduce a gravitational field which gives rise to mechanical forces induc- 
ing heat current flowJ-^ 

First, we, for simplicity, restrict the form of the BdG Hamiltonian to 



Aij is a constant Hermitian [i.e. (Al.) aiaV = (<'V/>, A \„ = <>V/),/w,vl matrix which preserves 
Ajj = Aj h and the first term includes the kinetic energy and the superconducting gap which 
includes second-order differential operators. C(r) is a Hermitian-matrix-valued function of 
r, and includes the periodic potentials, the superconducting gap, the Zeeman energy, and so 
on. Bj(r) is a matrix-valued function of r which is raised by gauge potentials, the spin-orbit 
interactions, or superconducting gap which includes first-order differential operators. Now 
we can get the equation: EiWO^i + «<9;fiJ(>)} = Z,{5,<>)/d ; + idiB t {r)) + [d u 5',(r)]}, where 
Bi(r) and 5',(r) are Hermitian matrices: B { (r) = 5,(r) + iB'j(r), and the last term of the RHS 
of the equation is noting but a Hermitian matrix, so it can be absorbed into C(r). Therefore, 
we take B t (r) Hermitian. Note that the Hamiltonians of the models mentioned above, the 
models of chiral p-wave (or rf-wave) superconductors and s-wave superconductors with the 
Rashba SOI, are expressed in the form eq. ©. From the the Hamiltonian eq. ©, we define 
the velocity operator : 



v,- is Hermitian and preserves PHS: v iasa ' S ' = +v* i _ as _ a , s >. (Note the sign "+" ) 

Next, we define the energy density operator and the energy current operator. By using the 



H BdG (r,d r ) = YjAijdidj + J^lBiirM + idtf]^)} + C(r). 



(5) 



ij i 



Vi = i [H BdG , n] . 



(6) 



method of integration by parts, the Hamiltonian ik can be rewritten as 

#= f d d rh(r), 
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where 



h(r) = - \-(dj{jj)\Ajjdjijj) + $\Biidi$) + + p(C$) 



2 

is the Hamiltonian density operator and Hermitian [i.e h'(r) = h(r)]. Here (diif/y(AijdjiJ/) is 
the abbreviation for 2Zijasa's'(di&as)H(Aij) asa ' S >d ji(/ a > s >}. A similar abbreviation is used for other 
terms. Therefore, in the presence of a gravitational field 0(r)J2E The Hamiltonian density 
operator and the Hamiltonian of the whole system are transformed into 

k^r) = (I + cf>(r))h(r), (7) 



where 



H<p(r, d r ) = ^ 



+ id i (l+ < P)B i (r)} + (l+ ( p)C(r)]. 

Note the scaling relation #<^_ = \H BdG . 

Now we define the energy current operator of the system with a gravitational field as 
follows, 

jE*i(r) = - -(1 + <f>(r)) {{v^{2H^) + h.c) 

- £ , 7 ^j(l+0(r)) 2 A,}], (8) 

where 

ol 

The last term of eq. © is indispensable for preserving the scaling law: 7 

jw(r) = (l + <Kry) 2 jEi<r), (9) 

where 

jEi(r) = j£^i(r)| 0=o = ^ ^(yi4rf(H BdG \f/) + h.c. - s ijk djk\ . 

We can check that the scaling law eq. © actually holds by a straightforward calculation with 
paying attention to Hf = (1 + <p)\H BdG - ^(d,-0)v,-. 

These Hamiltonian density and energy current operators eqs. © and © indeed satisfy 
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the continuity equation: 

—qT = [W.^J = -2j ■ do) 

i 1 

We present a brief proof of the continuity equation eq. (fTOl ) in the following. By noting the 
equations iff = -i [iff, < H^ = —1iH$ and v, = HAjjdj - 2B,, we can get 

^ =(1 + ^{-icVi^Cfl^) - (BAWiH+fo - i(Cfo\H+fo} 

+ h.c. . 

Note that the factor 2 of the RHS of the commutation relation shown above follows from the 
relation of operators \f/ as (r) = yL as (r) and the PHS of the Hamiltonian H^ asa > s> = -H* as _ a , s ,. 

On the other hand, by noting the equation <9,-(l + 0)v,- = 
M [Hp - ^Bi\ + cf>)dj - |(1 + 0)C], we get the equation: 

=i [{5,(1 + 0(r))v^} f (2H$) + ((1 + tir^ViWidilH^) 



+ /i.e. 

=(1 + ^(v^)^//^) + (Bidifo\H$) + i{C{ff)\H^)} 
+ h.c. . 

Therefore, we obtain the continuity equation eq. (TT0T >. 

Using the energy current operator eq. © satisfying the conservation law, we calculate the 
thermal Hall conductivity in the next section. 

Thermal Hall conductivity — In this section, we calculate the thermal Hall conductivity 
of superconductors using the procedure which was introduced by Qin, Niu and ShiP with a 
particular attention to the symmetry of the eigenfunction eq. © and compare the result with 
the case of normal metals and band insulators. 

The thermal Hall conductivity is given as follows: 

2M Z 

< = 4"'° + y#- (ii) 

The first term is given by the usual Kubo formula,^ and M Z E is the gravitomagnetic energy 
(heat) magnetization, which characterizes the circulation of the energy (heat) flow.^ 

From now on, we use a bra-ket notation: (f\0\g) means 
Rasa's' f d d rf* s (r)O asa , s ,(r,d r )g a , s <(r). 

Now we introduce new operators, which are the Fourier transforms of the field operators, 
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and expand them in the operators fib, and ^ kn : 

tLg = I d d re~ K ~ q)r h(r) = ^ fy\J/k>n>h~ q knk>n' , 

J knk'ri 

jEqi = I d d re~ lqr j E i(r) = ^ $l„fa>n> jEqxknk'n' 
^ hi k'n' 



where 



i - 1 / f i H BdG^' q ' + e iqr H BdG 2 

n-qknk'n' = -j\Jkn\ ~ ^ k ' n '' > 

(as q 0), (12) 

(H BdG fkn\e' k "-\Vjf kw ) + (Vjf kn \e~ k "-\H BdG f kw ) 



JEqjknk'n' — ^ 



8/ 

i 

-(d iVi f kn \e~ iqr \vjf k ^) - (vifknle-^ldiVjf^)}] , (13) 

which can be obtained by noting the two identities: Ajjdie iqr dj = (Aijdidje lqr + e lqr Ajjdidj)/2 + 
0(q 2 ) and e^Bjidi + idfifi^ = {e~ iq, \Biidi + idfa) + (B t id t + idiBde^}/!. Note that the 
coefficients for the expansion, h- qknk ' n > and jEqjknk'n', preserve the symmetries 

h-qknk'n' h—q —k'-n' —k—n 5 

jEqjknk'n' j E q j — f —n' —k—n> (14) 

which follow form H BdGasa , s , = -H* BdG _ as _ a , s , and v iasa , s > = v?_ fl „_ flV , and eq. ®: f kn as(r) = 
f*t (r). 

The first term of eq. (TTTT > is give by 

1 r°° 

4"*° = -^J o dte~ Qt (J Ey ;J Ex (t)), (15) 

where 7 £( - = j£ ? /| ?=0 , = l/B dA(a(-iA)b), d(t) = e^'de" 1 ^ 1 ' , and V is the volume of 
the system. By the formula for the four-point correlation function eq. © and the symmetry 
of the coefficient eq. (fT4l) . we get 

Kubo _ 1 f(Ekn) ~ f(Ek'n') 

K *y ~ vt 2 U 



i(E kn - £W) 

knk'rr 
(kn)^(k'n') 



2 



X J Eyknk'n'(.J Exk'n' kn JEx—k—n—k'—n') 



~ x rrrt'} / , . , nt J Ev kn k' n' J Exk' n' km 

VT 2 i(E k „ - Ev„,) 2 



kn k'n* 
(kn)^(k'n') 



,r, j-i nt Evknk'n' Ex k'n' km 

yt^kn ~ t^k'n') 



8/DU 
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where JEiknk'n' = jEqjknVn'\ =0 - The factor 2 in front of eq. (1161 ) is a result of the PHS. 
Moreover, by calculating in the manner similar to ref. [15] we obtain 



K 



Kubo 



—y 

4TV ^ 



Im 



kn 



du 



kn 



dk x 



(H B dGk + Ekn) 2 



du 



kn 



dh 



f(Ekn), 



(17) 



where H BdGk = e ,kr H BdG e 



„ikr 



Next we calculate the gravitational magnetization M Z E . It is the solution of the differential 
equation 



2M% - T 

E dT 



dM% J3 ( d , , d , , 

~-2i\W x { q,jEqy) ~d^ { q,jEqx) 



with a boundary condition Um T ^oT -jf- = 0. In order to evaluate it, we also carry out a 
calculation similar to ref. [T5] with paying attention to the last extra term of eq. © and the 
symmetries eq. (fT4l) . We get 

du kn ' 



kn 





ldu kn 






\ dk x 



du 



(HBdGk + E kn ) 

du kn 



3ky 



f(E kn ) 



kn 



dk x 



dky 



+4Imi 



1 du 



ki, 



dk x 



du kn \ r 

dky I X 



f(Ekn) 

Ekn 

xf(x)dx 



(18) 



Note that the expressions eqs. (fTTT) and (fT8l are one half of those in the case of normal 
metals or insulators [cf: eq. (23) of ref. |7]and eq. (S97) of ref. Q3). This is due to the factors 
1/2 of eqs. O and O [cf: eqs. (S75) and (S76) of ref. E] and 2 of eq. flU) [cf: eq. (S60) of 
ref. [151, which are caused by the last extra term of eq. @ and the relations eq. (fT4T) associated 
with the PHS. 

In the end, we obtain an expression for the thermal Hall coefficient: 



4 = -^/^I>fe 



kn 

Ekn<E 



du h 



kn 



dh. 



f'(E). 



(19) 



Moreover, by using the Sommerfeld expansion, we obtain an expression in the low- 
temperature limit: 

where C\(E) is the TKNN number, which is an integer when the energy E lies in the en- 
ergy gap, 132 and it is given by, C X (E) = £ n f v Im «*) ®( E ~ Ekn ^ where 0(jc) is the 
Heaviside step function. 

It is notable that the quantization value i^Ci(O) is exactly one half of the value of the 



9M 



J. Phys. Soc. Jpn. 



LETTERS 



Chern insulator (or the IQHE state). In the case of spineless chiral p-wave superconductors, 
the TKNN number Ci(0) is equal to ±lj^ and thus, this result is in agreement with the result 
obtained from the Ising CFT with central charge c = 1/2 for the edge state 

Summary — We have demonstrated that the thermal Hall conductivity of 2D TSCs with 
broken time reversal symmetry is quantized: K xy = -y^f"' where C\ is the TKNN number of 
the BdG Hamiltonian. Our approach, solely, relies on bulk calculations, without referring to 
the Majorana edge theory. This value is one half of the thermal Hall conductivity in 
the case of the Chern insulator (or the IQHE state). In the derivation of this result, the PHS, 
which BdG Hamiltonians generally possess, plays an important role. Our result is in perfect 
agreement with that obtained from the Ising CFT which describes the edge state of 2D TSCs. 

We are thankful to Ken Shiozaki for helpful discussions. This work is supported by the 
Grant-in- Aids for Scientific Research from MEXT of Japan (Grants No. 23102714 and No. 
23540406). 



lOffi] 



J. Phys. Soc. Jpn. 



LETTERS 



References 

1) D. J. Thouless, M. Kohmoto, M. P. Nightingale, and M. den Nijs: Phys. Rev. Lett. 49 
(1982) 405. 

2) M. Kohmoto: Ann. Phys. (N.Y.) 160 (1985) 343. 

3) A. P. Schnyder, S. Ryu, A. Furusaki, and A. W. W. Ludwig: Phys. Rev. B 78 (2008) 
195125. 

4) G. E. Volovik: The Universe in a Helium Droplet (Clarendon Press, Oxford, 2003), 
p.136. 

5) N. Read and D. Green: Phys. Rev. B 61 (2000) 10267. 

6) K. Nomura, S. Ryu, A. Furusaki, and N. Nagaosa: Phys. Rev. Lett. 108 (2012) 026802. 

7) T. Qin, Q. Niu, and J. Shi: Phys. Rev. Lett. 107 (201 1) 236601. 

8) In ref. [6l it is shown that the thermal Hall conductivity of spinfull Majorana fermions 
on the surface of a 3D TSC is half-qunaized and associated with the TKNN number. 
Moreover, this model is related to that of spinless p + ip superconductors via a unitary 
transformation: c = (fa + i^j,)/ V2, c 1 " = (fa - iiff])/ V2. 

9) M. Sato, Y. Takahashi, and S. Fujimoto: Phys. Rev. Lett. 103 (2009) 020401. 

10) M. Sato and S. Fujimoto: Phys. Rev. B 79 (2009) 094504. 

11) M. Sato, Y. Takahashi, and S. Fujimoto: Phys. Rev. B 82 (2010) 134521. 

12) A. Altland and M. R. Zirnbauer: Phys. Rev. B 55 (1997) 1142. 

13) J. M. Luttinger: Phys. Rev. 135 (1964) A1505. 

14) The chemical potential of the Bogoliubov quasiparticles is measured from the Fermi 
level, and set equal to 0. Thus, the energy magnetization is equivalent to the heat mag- 
netization. 

15) Supplemental Material of ref. [7] 



